Acoustic vibrations of a silica-embedded gold nanoparticle: elastic anisotropy 
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Spherical gold nanoclusters are grown by 1.8 MeV Au ++ ion implantation into an amorphous silica 
matrix and subsequent air annealing at 873 K. Ultraviolet and visible light absorption confirms the 
presence of a dipolar surface plasmon peak at the expected location for a gold sphere in silica. 
Grazing incidence X-ray diffraction shows peaks corresponding to fee gold with lattice spacing close 
to that of bulk gold, as well as allowing estimation of cluster size using Scherrer's formula. Low 
frequency Raman scattering reveals a relatively narrow peak, suggesting a narrow distribution of 
nanocrystal diameters. Acoustic phonon frequencies corresponding to the spheroidal quadrupolar 
vibrations of a continuum sphere with the anisotropic elasticity of gold are calculated using a novel 
method of molecular dynamics and extrapolation to the continuum limit using relatively small 
numbers of point masses. The study confirms high energy ion implantation as a method capable of 
producing gold nanocrystals with spherical shape and well controlled size. 



I. INTRODUCTION 

A free elastic sphere has vibrational modes whose low- 
est frequency is on the order of the speed of sound divided 
by the sphere radius. This classic problem of elastic me- 
chanics was solved by Lamb in 1883i for the case of a 
homogeneous sphere with isotropic elasticity. 

Ninety years later, this solution found application in 
explaining anomalies of the specific heat of powders of 
spherical lead particles 2 to 4 nm in size.- Direct obser- 
vations of vibrational modes of mixed spinel (MgC^Ov 
MgAi204) microcrystallites 22 to 35 nm in diametes^ be- 
came possible by shining an intense laser beam onto the 
sample and carefully analyzing the scattered light to de- 
tect scattered frequencies slightly above and slightly be- 
low the laser frequency. The frequency shift of the light 
is equal to a vibrational frequency of the spinel micro- 
crystallite. Since then, Raman and Brillouin scattering 
have remained powerful tools in the characterization of 
"nanoparticles" - particles with diameters on the order 
of a few to a few dozen nanometers. 

The nanoparticles in these experiments are often not 
free, but rather embedded in solid materials such as 
glass. The Lamb solution does not directly apply since 
it assumes zero traction force boundary conditions of 
the sphere (i.e. a free surface). Initial attempts^ to 
generalize the calculation to an isotropic elastic contin- 
uum sphere embedded in an infinite clastic matrix are in 
error. 6 But an earlier calculation 7 in the geophysics liter- 
ature had correctly solved the problem and showed that 



embedded sphere vibrations can be acurately approxi- 
mated with complex-valued pseudo-frequencies w where 
Re(w) is the center frequency of the Raman peak and 
2Im(w) is the full width at half maximum (FWHM). 
Correct application of this solution to nanoparticles was 
made only recentlyi&Siifli 

Vibrations of an embedded sphere are described in 
terms of a displacement field u(f , t) which equals r ' — R 
where R is the equilibrium position of a material point 
and r is its displaced position. Due to the spherical sym- 
metry, it is convenient to choose the origin at the center 
of the sphere and use spherical coordinates r, and 4>. 
The time dependence is exp(— iujt) where a; is a complex 
number. In general, a vector field u is the sum of a zero- 
curl field V$ and a zero divergence field VxH where 
$(r) is a scalar field and S(r) is a vector field. 

In the situation of an isotropic, homogeneous elastic 
medium with transverse speed of sound vt and longi- 
tudinal speed of sound v^, normal modes are found by 
solving two uncoupled equations 

(v 2 + fc!)$ = o (i) 

and 

(V 2 + 4)S - (2) 

where kr, = uj/vt, and kx — uj/vt- 

For convenience, the original Lamb solution will be 
called the Free Sphere Model (FSM). The isotropic 
sphere embedded in an infinite matrix will be called the 
Complex Frequency Model (CFM). 
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Whether for FSM or CFM, there are some solutions 
of the equation of motion with V • u = 0. Lamb called 
these the "first class" and the modern word is "torsional" 
(TOR). Apart from torsional modes, an important spe- 
cial case are radial or "breathing" modes, in which u is 
purely along the r-direction and depends only on r and 
not on 9 or <fi. In this case u is of the form V$. The re- 
maining modes of what Lamb called the "second class" , 
now called "spheroidal" (SPH), have a more complicated 
form: 

m = V$ + VxVx (rip) (3) 

where $ and ip are scalar fields. Apart from the breathing 
modes, spheroidal modes are neither zero curl nor zero 
divergence. 

It can be said that torsional modes are purely trans- 
verse in nature. As a result their frequencies depend only 
on vt ■ The breathing modes are purely longitudinal, but 
mode frequencies depend on both vt and vl- The re- 
maining spheroidal modes are a mixture of motions, but 
in most cases their frequencies depend primarily on vt 
and only weakly onwi, so they are predominantly trans- 
verse in character. For a generic material with Poisson 
ratio |, 80% of all modes, counting both TOR and SPH, 
are more transverse, while 20% are more longitudinal. 

To systematically label FSM or CFM modes, it is con- 
venient to use index q e {TOR, SPH}. In addition, modes 
have an angular momentum number £ £ {0,1,2...} and 
a z-angular momentum m where — £ < m < I. For 
brevity we consider only the m = case below. But it 
is important to keep in mind that FSM and CFM modes 
have degeneracy 21+1. There is also a solution index 
n e {0,1,2...}. The scalar field $ satisfying Eq. JIJ has 
the form 

§(r,0)=j t (k L r)Pt(coB6) (4) 

where ji are spherical Bessel functions and Pe are Leg- 
endre polynomials. The scalar field ip in Eq. © has the 
form 

iP(r,e)=jt(k T r)P e (cos9). (5) 

A given normal mode, whether FSM or CFM, will be 
denoted (q,£,m,n) and its frequency is or w^ M . 

The first applications of the correct CFM solution 
for an isotropic continuum sphere embedded in a ma- 
trix were restricted to spheroidal modes Group the- 
oretical arguments^ show that in the absence of any 
kind of anisotropy or non-sphericity, and in the limit 
that the nanoparticle is much smaller than the wave- 
length of the laser light, only the modes (SPH/=0) and 
(SPH,i?=2) can lead to Raman scattering. Correct calcu- 
lations of the CFM torsional and spheroidal £ ^ mode 
frequencies were reported lateral Even so, the usage of 
complex- valued frequencies of normal modes is problem- 
atic. There is no corresponding displacement field. 

Clarification of the significance of CFM calculations 
was obtained through an exhaustive calculation of all 



of the real-valued frequencies of an isotropic spherical 
nanoparticle embedded in a large, but finite matrix with 
a free spherical outer surface of radius R m & This is called 
the Core-Shell Model (CSM). Taking the limit R m -> oo, 
the number of CSM modes in a given range of frequency 
becomes infinite, and has the Debye density of states. 
However, assuming uniform random phase excitation of 
the CSM modes, the mean squared displacement within 
the nanoparticle exhibits peaks whose positions agree 
closely with the real part of CFM modes, and whose 
FWHM agree closely with twice the imaginary parts of 
CFM modes. Thus, CFM calculations have a solid basis 
and utility, even if the corresponding CFM modes must 
be treated with caution. 

Both FSM and CFM have limitations of applicability 
to real systems. (1) Nanoparticles are not truly spherical 
in shape, although the dynamics of most elaboration pro- 
cesses make them approximately that way and transmis- 
sion electron microscopy (TEM) images confirm this. (2) 
For small nanoparticles (smaller than those in the present 
study), the wavelength of even the low- lying modes is not 
long compared to interatomic spacing, so that we are not 
close to the center of the Brillouin zone and the acoustic 
limit. (3) There will be deviations of the properties of the 
nanoparticle material from the same material in the bulk 
due to surface effects. (4) In a given sample, the radii 
of nanoparticles will have a distribution, rather than a 
single value. (5) The interface between the nanoparticle 
and the matrix is not always correctly described in terms 
of continuity of the displacement field and components 
of the stress tensor^ (6) The speed of sound within the 
nanoparticle depends on the direction of propagation if 
the elasticity of the material is anisotropic. Only the last 
issue is addressed in this paper. 



II. ELASTIC ANISOTROPY 

Since these words are frequently misused, we recall 
that "isotropic" means "unchanging under rotation" 
while "homogeneous" means "unchanging under trans- 
lation" . Although bulk polycrystalline gold is elastically 
isotropic, a single crystal of gold is not. In fact, the no- 
ble metals (Au, Ag, Cu) are remarkably anisotropic. De- 
pending on the direction of propagation in the crystal, a 
transverse acoustic (low frequency) sound wave can travel 
at speeds varying from 867 m/s to 1482 m/s. The longi- 
tudinal speed of sound varies less, ranging from 3147 m/s 
to 3440 m/s. But since all FSM and CFM modes depend 
on vt, it must be expected that the strong transverse 
anisotropy will play an important role. 

Application of FSM to elastically anisotropic crystals 
(in fact, no crystalline material is elastically isotropic) 
goes back many years. Since the theory asked for vt 
and vl, it was common to simply use the speeds of 
sound in the bulk polycrystalline material, sometimes 
with apologyi^ but usually without. One approach is 
to make use of the speeds of sound along the principal 
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crystal directions, (e.g. <100>, <110> and <111> in 
siliconi4) so as to get several different estimates of the 
frequencies. 

The elastic mechanical problem of calculating the vi- 
brational modes of a free sphere with anisotropic elas- 
ticity can be set up analytically for certain special kinds 
of anisotropy^ but not for cubic crystalline elasticity 
such as that of fee gold. Numerical techniques for find- 
ing eigenfrequencies of an object with arbitrary elasticity 
and a variety of special shapes have long been knownii§ii£ 

Our numerical method used here varies only slightly 
from a previous one^S Its features are: (1) arbitrary 
cubic crystalline elasticity including arbitrarily high 
anisotropy (unlike the earlier method 10 which numeri- 
cally blows up when faced with the anisotropy of a no- 
ble metal, but worked for silicon); (2) ease of handling 
a nanoparticle of arbitrary shape; (3) ready extendabil- 
ity to inhomogeneous density and elasticity; (4) gener- 
alizability (with simple modification) to hexagonal and 
other crystal symmetries; (4) self-diagnosis of numerical 
convergence; (5) rapid convergence; (6) conceptual sim- 
plicity. 

Our approach can be described as a Finite Difference 
Time Domain (FDTD) calculation, and also as molecu- 
lar dynamics (MD), the former phrase being common in 
the engineering literature. The continuum sphere of ra- 
dius R p and density p p is approximated by N equal point 
masses arranged in a simple cubic (SC) lattice. The point 
masses do not represent atoms. For calculations reported 
here, N was as high as 8000. The use of an SC lattice is 
associated with the cubic symmetry of the elasticity but 
does not preclude simulations of objects of crystalline 
material with fee, bee or diamond structures, as well as 
isotropic solids such as glass. The method as detailed 
here applies only to the homogeneous cubic case. The 
method is unrelated to calculations which model the mi- 
croscopic inter-atomic forces such as the Tersoff potential 
for silicon or the Keating model. Rather, it is a contin- 
uum model of the material, involving only the bulk elastic 
constants and the density. 



III. MOLECULAR DYNAMICS METHOD 

Our method is a molecular dynamics calculation. 
While this approach can readily be extended to the more 
general case of an object of arbitrary shape with inho- 
mogeneous anisotropic elasticity, we choose for brevity 
only to present the case of a spherical object with cubic 
crystalline elasticity. A simple cubic lattice with lattice 
parameter a of N identical point particles of mass pa 3 in- 
teracting through "springs" has particle positions fj inte- 
grated in time and Fourier transformed in order to obtain 
mode frequencies. Each particle has six first neighbors, 
twelve second neighbours, eight third neighbors and so 
on. Every pair of first neighbors are coupled by po- 
tential energy \k sp i(\\ri—rj\\—a) 2 . Second neighbor pairs 
are coupled by \k sp 2(\\fi — r } -\\ — \/2a) 2 . Third neigh- 



bor pairs are coupled by \k sp3 ( 



Fourth neighbor pairs (i, j) are coupled by |fc sp 4(||ri — 
fj\\ — V4a) 2 . Fifth neighbor pairs (i, j) are coupled by 
2" ^sp5 ( 1 1 — — \/5a) 2 • Pairs are not double counted. 
Cubic octets are coupled by (ks, p t/^a 2 )(D — 6a 2 ) 2 where 

D = J2j=i \\Rcm - r-jW 2 and R cm = ± J2j=i 

No microscopic physical interpretation is offered for 
the ks p t term. Its purpose it to allow reproduction of 
the bulk elastic constants rather than to represent the 
potentials due to atomic bonding among atoms. It is 
possible to approximate interatomic interactions among 
atoms in certain special kinds of crystals using model 
potentialsiiSiiSiSS Our approach is not limited to certain 
kinds of crystals, but rather to continuum materials with 
completely general cubic elasticity. For this purpose, it 
is numerically more efficient to use the potentials that we 
employ. 

Although we have introduced six force constants, (k sp i, 
k sp 2, k sp3 , k sp4 , k sp5 , and k Spt ) only k spl , k sp3 , and k 8pt 
are used for the numerical calculations of this paper. A 
previous paper— used k sp i, k sp 2, and fcg pt , but was lim- 
ited to materials of lower elastic anisotropy as will be 
explained below. It is advantageous to use the earlier 
method for less isotropic materials since numerical con- 
vergence is faster as N increases. 

The next step is to show how the force constants are 
related to the bulk elastic constants of the material of the 
nanosphere. A cubic crystal (in particular, any material 
with bec, fee, diamond or zincblcnde crystal structure) 
has three elastic constants: Cn, Cyi, and C44. For an 
isotropic material, C12 = C\\ — 2C44. The degree of 
elastic anisotropy is quantified by the Zener anisotropy 
ratio, Az, where 



A z = 



2C, 



44 



Cn — C12 



(6) 



Elastic isotropy corresponds to Az = 1, and anisotropic 
materials can have Az above or below 1. The Zener 
anisotropy of gold is 2.92. Silver and copper have com- 
parable anisotropies. 

The bulk modulus of an arbitrary cubic elastic mate- 
rial is B = (Cn + 2Ci2)/3. Let V denote volume. A 
small isotropic dilatation S = dV/V stretches all neigh- 
bor distances proportionally and leads to stored energy 
^VB5 2 , so that we obtain 

a(Cii+2Ci2) = k spl +4:k sp 2+4:k sp3 +4:k sp 4+20k sp5 +24:ks p t 

(7) 

Considering that a longitudinal plane wave moves 
along the x-axis at speed ^JUuJp, we find 



aCn = k sp i + 2k sp2 - 



-4k sp 4, + 13M sp5 +8k Spt (8) 



Finally, considering that a transverse plane wave moves 
along the x-axis at speed ^JC^/ p, we find 



aC*44 



%j2 



3.2fc 



Sp5 



(9) 
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By subtracting Eq. [5] from Eq. [7| we obtain 
4 

aCx2 = k sp 2 + -^k sp3 + 3.2k sp5 + 8k 8pt (10) 

Note that if k% pt — 0, Cyx and C44 are exactly the 
same. This has been numerically verified up to tenth 
neighbours, ft is coincidentally the case that Cn and C44 
are almost the same for NaCl, so the elastic properties of 
NaCl could be modelled without the kg, pt term. In addi- 
tion, C12 and C44 are equal for isotropic materials with 
Poisson ratio 4, such as polycrystalline Cerium, once 
again allowing modelling with two-body forces alone. 

For other materials, it is essential to include an interac- 
tion apart from two-particle forces. This is the reason for 
including the ks p t potential term. Eight masses must be 
involved in order to have a term that in itself respects 
cubic symmetry. However, it would also be workable 
though not obviously advantageous to have used a three 
body interaction term. 

Three interaction terms, including k% pt , must be used 
to exactly fit the three elastic parameters Cn, C12, and 
C44 . If, as in earlier workifl, only k sp \ , k sp 2 , and ks P t are 
used, then 

k sp i = a(Cu - (C12 + C44)), (11) 

Note that materials with Zener anisotropy above 2 will 
have negative k sp \ . This leads to instability of molecular 
dynamics simulations of the lattice. 

Another option is to use only k sp \, k sp3l and kg pt . In- 
version allows the three force constants to be obtained 
for a general cubically elastic material: 

k sp i = a(Cu - C 12 ), (12) 



k S p3 — —t-Caa (13) 

and 

k Sp t = \a{C 12 - C 44 ). (14) 
o 

Note that k sp \ and k sp s are never negative. It can be 
the case that k 8pt is negative for some materials, but in 
practice this does not lead to instability of the molecular 
dynamics simulations. As a result, any cubic material 
can be simulated. 

The frequency dependence of a given mode is f(N). 
A sequence of simulations with N varying from a few 
hundred to 8000 is used to extrapolate /(oo), which is the 
continuum limit. Accurate results are obtainable with 
just a few thousand. 

The correctness and degree of convergence of our MD 
method was checked by comparing its results to those 
of the Lamb solution for isotropic elastic materials of 
various Poisson ratios. The results of the two methods 
agree closely. 



TABLE I: Frequencies of selected vibrational modes of an 
elastic continuum sphere with the elastic constants of single 
crystal gold of radius R are shown. The value given is uiR (in 
m/s) where to is the frequency in radians per second and R is 
the nanosphere radius in metres. 
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IV. MOLECULAR DYNAMICS RESULTS 

A sequence of simulations with increasing N is per- 
formed in order to extrapolate the limit of infinite N . 
The purpose of the calculation is to solve for the vibra- 
tional frequencies of a continuum elastic sphere, and does 
not take into account the crystal structure of the mate- 
rial or dispersion (frequency dependence of sound speed) 
of phonons. 

The time duration of the simulation is determined by 
the frequency resolution that is desired. The initial ve- 
locity distribution corresponds to an FSM mode. The 
z-axis of the desired FSM mode is randomly oriented for 
each run relative to the axis of the cubic lattice. 

A Fourier transform is carried out of the motion of 
several randomly selected mass points along various axes. 
A plot is then made of the power spectrum versus 1/N. 
This permits extrapolation of the frequency of each mode 
to infinite N. The extrapolation is done visually. 

The elastic constants usedSi for monocrystalline fee 
gold at 300 K were Cn = 191 GPa, C 12 = 162 GPa, 
C 44 = 42.4 GPa, and p = 19.283 g/cc. 

Tablc[|]shows results of the molecular dynamics simula- 
tions. Frequencies are presented in terms of the ratio loR. 
For an isotropic material, it is beneficial to use dimen- 
sionless frequencies such as 77 = loR/vt and £ = loR/vl- 
In this problem, the anisotropy of the elasticity means 
that there is no convenient sound speed scale to form a 
dimensionless frequency. However, expressing results in 
terms of u>R exploits the scale invariance of the problem 
as long as continuum elasticity is assumed. 

The frequencies in Tab.[I]are classified according to the 
symmetry of the initial disturbance used to excite them. 
In each case, velocity fields of spheroidal FSM modes of 
angular momentum £ are used as the initial disturbance. 
The z-axis of the velocity field is randomized with re- 
spect to the axes of the lattice for each run. The mode 
degeneracies in Tab. [I] were found by gradually distort- 
ing the shape of the sphere so as to break all symmetry, 
and noting into how many frequencies each mode splits. 
The £=0 spheroidal mode frequency is calculated, but 
the Raman scattering strength is typically much smaller 
than for the 1=2 spheroidal mode. 

It would be expected on the basic of group theoretical 
considerations that a 5-fold degenerate (SPH, 2) mode of 
an isotropic sphere would be broken into two different 
levels by lowering the elastic symmetry to that of a cu- 
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TABLE II: FSM and CFM frequencies for a directionally 
averaged (elastically isotropic) sphere of gold. For the com- 
plex CFM frequencies, the gold sphere is embedded in a silica 
matrix. The value given is uoR where uo is the frequency in 
radians per second and R is the sphere radius in metres. 
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FIG. 1: Perturbed frequencies of the initially fivefold de- 
generate quadrupolar spheroidal modes of a free sphere are 
shown as the parameter x is varied, x — corresponds to 
directionally averaged isotropic gold while x = 1 corresponds 
to true anisotropic gold. The upper and lower branches are 
twofold and threefold degenerate respectively. 



bic crystal. The corresponding representations are E g 
with two dimensions and Ti g with three dimensions. 22 
To illustrate this, Fig. ^ shows how the frequencies vary- 
as the amount of anisotropy is continuously varied. Let 
Ciia, Ci2a, and C^ a be the elastic constants of direction- 
ally averaged gold. These are based on the directionally 
averaged sound speeds of 3330 m/s and 1250 m/s. In 
particular, C lla = (3330 m/s) (3330 m/s) (19293 kg/m 3 ) 
and C 44 a = (1250 m/s) (1250 m/s) (19293 kg/m 3 ). Also, 
since isotropy is assumed, C\2 a = G\\ a — 1Cm a . The elas- 
tic constants of true monocrystalline bulk gold are Cm,, 
Ci2h, and C44h as given earlier in this section. The fre- 
quencies in Fig. ^ are done with elastic constants CV, = 
xCij a + (l — x)Cijb- Note how the isotropic (SPH, 2) mode 
at loR — 3320 m/s continuously splits into two modes. 



V. MATRIX EFFECT 

The molecular dynamics calculation results shown in 
Tab. [I] assume a gold sphere with a free surface. However 
the gold nanospheres reported in the experimental results 
of this paper are embedded in amorphous silica. As a 
result, the outer surface of the gold spheres is not free. 
The silica will affect the vibrational frequencies. 

The complex frequency model (CFM)£i2i±!i calculates 
the complex-valued vibrational frequencies of a sphere 
embedded in an infinite matrix. Both the sphere and 
the matrix are assumed to be continuous, homogeneous 
and isotropic. Since vibrational modes of the gold sphere 
are damped by the propagation of sound waves into the 
matrix, the modes have complex valued frequencies. The 
imaginary part of the frequency is related to the damping 
time of the mode. 

To treat gold as an isotropic material requires an ap- 



proximation to be made about its elastic properties. We 
use directionally averaged speeds of sound. The longitu- 
dinal and transverse speeds are separately averaged. The 
results are 1250 m/s and 3330 m/s. 

Amorphous silica is an isotropic material. We assume 
speeds of sound of 3760 m/s and 5950 m/s and a density 
of2.20g/cc. 

Results of CFM calculations are shown in Tab. El FSM 
calculations based on Lamb's 1882 solution for a free 
sphere are also given. 

We can either model the anisotropic elasticity of gold 
using molecular dynamics calculations for a free sphere or 
else model the effects of the silica matrix using CFM for 
an isotropic sphere. We have no way to simultaneously 
model both anisotropy and matrix effects. 

In Tab. El note that the damping of the modes is rel- 
atively small, based on the low imaginary parts of the 
CFM frequencies. However, the (SPH,^=2) mode is sign- 
ficantly upshifted by the silica matrix. This is caused by 
the action of the matrix to provide a spring-like reaction 
force acting at the outer surface of the gold sphered As 
a result, the MD calculations which ignore the matrix 
effect will require some correction. 



VI. TANDETRON ACCELERATOR 

Gold ion implanation was done using a 1.7 MV Tande- 
tron accelerator (Model 4117HC) from High Voltage En- 
gineering Europa, Amersfoort, the Netherlands. This ac- 
celerator has the "tandem" configuration where negative 
ions generated by a cesium sputter gold ion source are 
first accelerated from ground to the high voltage termi- 
nal and there the negative gold ions are converted to 
positive gold ions, while passing through a stripper canal 
filled with nitrogen gas. The same high voltage again 
accelerates the positive Au ++ ions to ground potential. 
In this way, Au ++ ions could be generated with energies 
up to 5.1 MeV with the Tandetron. 



VII. SAMPLE PREPARATION 

Using the Tandetron ion accelerator, Au clusters were 
grown by direct Au-implantation on fused quartz (amor- 
phous silica) substrates. Implantation was performed us- 
ing 1.8 MeV Au ++ at lxl0~ 5 Pa with a beam current 
of wl2 mA m -2 to avoid substantial beam heating in 



6 



the ion fluence range of 2xl0 20 - lxlO 21 m~ 2 . Dur- 
ing implantation the substrate was at room temperature. 
The diameter of the ion beam is 1 cm, so that the to- 
tal beam current is approximately 1/iA. At this energy, 
the average penetration depth of the Au ++ in the silica 
is 445.6 nm before they come to rest, with a straggling 
(distribution width) of 55.4 nmiSi The advantage of using 
such a high beam energy is that structure can be created 
at a much greater depth below the surface of the silica 
substrate. High temperature air annealing at 873 K for 
1 hour was performed to promote growth of the clusters 
as gold atoms become mobile in the silica substrate and 
are able to conglomerate into bulk gold metal as opposed 
to isolated gold atoms. 

Our study emphasizes acoustic phonons of the im- 
planted nanospheres and therefore benefits from a nar- 
row size distribution. The 1.8 MeV beam energy is not 
thought to play a role apart from controlling the implan- 
tation depth of the ions. Ion implantation often results in 
wide size distributionsi 25 i 26 Implantation temperature 2 ^ 
and post-implantation annealing 2 - are known to play an 
important role in controlling average particle size and size 
distribution. The current density of the ion beam affects 
average size and size distribution when implanting silica 
with Ag2£ and Cu^i and is expected from simulations 2 ! 
to play a similar role for gold. Laser pulse annealing^! is 
also known to be able to narrow the size distribution. 



VIII. SAMPLE CHARACTERIZATION 

Optical absorption spectra were recorded at room tem- 
perature in the range of 200 - 1100 nm using a Hewlett- 
Packard diode-array UV-visible (UV-Vis) spectropho- 
tometer (Model 8453) with necessary signal correction for 
the substrate. Figure [21 shows broad peak structure cor- 
responding to the frequency of the dipolar surface plas- 
mon resonance of a gold sphere in a silica matrix. This 
is evidence of the shape of the nanoparticles being ap- 
proximately spherical. The surface plasmon resonance of 
an ellipsoidal sphere would split into higher and lower 
frequencies. 

Grazing incident X-ray diffraction (GIXRD) was done 
using Cu Ka X-rays (wavelength = 0.15406 nm), as 
shown in Fig. 01 Both the (111) and (200) peaks were 
observed, with (111) being the stronger feature. The lo- 
cation of the (111) peak closely matches that expected 
based on a fee unit cell size of 0.40786 nm for gold at 
room temperature. This was done in 0.05 degree steps 
with 25 s accumulation time per data point. This is an 
£1-26 scan, fi is 0.8 degrees. In other words, the incident 
X-rays enter the sample at 89.2 degrees from normal in- 
cidence. 

X-ray diffraction is an indication of the crystalline do- 
main size rather than the size of the nanosphere. How- 
ever since annealing at 873 K was done, it is expected 
that the nanosphere would consist of a single crystal do- 
main, and the size from GIXRD is interpreted as the size 
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FIG. 2: Light absorbance as a function of photon energy 
is shown for post-annealed Au nanoclusters in silica matrix 
grown at fluence lxlO 21 m -2 . Structure around 2.4 eV is 
attributed to the dipolar surface plasmon resonance of a gold 
sphere embedded in silica. 
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FIG. 3: Grazing incidence X-ray diffraction data is shown 
for gold nanospheres embedded in silica, using copper alpha 
X-rays. Two expected peak positions are indicated based on 
the fee structure and 0.40786 nm lattice constant of bulk gold. 
The indicated averaged diameter estimated is based on Scher- 
rer's formula applied to the (111) peak. 



of the nanosphere. 

Low-frequency Raman scattering studies were done at 
room temperature in the back scattering geometry, us- 
ing vertically polarized 488 nm line of an argon ion laser 
(Coherent, USA) with 200 mW power. Unpolarized scat- 
tered light from the sample was dispersed using a dou- 
ble monochromator (Spex, model 14018) with instrument 
resolution of 1.4 cm" 1 and detected using a cooled pho- 
tomultiplier tube (FW ITT 130) operated in the photon 
counting mode. See Fig. 0] 

A strong and relatively narrow Raman peak is seen 
in Fig. 0] with its peak at 22 cm" 1 . The full width at 
half maximum (FWHM) of this peak is approximately 



7 



CO 




10 20 30 40 50 60 
Raman Shift (cm -1 ) 



FIG. 4: Room temperature low-frequency Raman spectra for 
post-annealed Au nanoclusters in a silica matrix. The 4 cm - 
wide peak at 22 cm -1 is attributed to thermal spheroidal 
quadrupolar vibrations of the gold sphere. The broadening of 
the peak may be primarily due to damping of the vibration 
due to mechanical contact with the silica matrix. 

4 cm -1 . This is relatively narrow compared to Raman 
peaks seen in other systems of embedded nanoparticles. 
The center of the Raman peak corresponds to an angular 
frequency ui = 4.15xl0 12 rad/s. 

IX. DISCUSSION 

The low frequency Raman spectrum of our sample in 
Fig. 0] showed a prominent peak with an angular fre- 
quency of 4.15xl0 12 rad/s. The average radius of the 
particles, as estimated from GIXRD in Fig. fusing Scher- 
rer's formula was 0.90 nm. The resulting experimental 
value of ojR is 3732 m/s. 

The existence of a vibrational mode does not mean 
that it will appear in a Raman spectrum. For an isotropic 
sphere, only the (SPH,0) and (SPH,2) modes can be Ra- 
man active based on group theoretical considerations^ 
and in practice the intensity of the (SPH,0) is much 
smaller. More detailed analysis, not presented here, is 
necessary to predict the relative Raman activity of the 
two different frequencies into which the (SPH,2) mode 
has been split by elastic anisotropy^ 2 ^ 

It is expected that one of the two modes with ujR = 
2342 m/s and 3779 m/s as shown in Tab. Q] will be the 
dominant feature in the Raman spectrum, but these val- 
ues from MD omit the matrix effect. To estimate the 
effect of the matrix, we rely on comparison of FSM and 
CFM results in Tab. |Hj which shows that the matrix ef- 
fect increases the frequency of the (SPH,£=2) mode by a 
factor of 1.28. Applying this to the MD results, the two 
expected locations of Raman peaks are u>R — 3006 m/s 
and 4837 m/s. One is 20% lower than the location of the 



experimental Raman peak while the other is 30% higher. 

Based on the cubic symmetry, for modes with Ti g sym- 
metry, the Raman tensor is off-diagonal while for modes 
with Eg symmetry, the Raman tensor is diagonal^ In 
other words, Ti g modes will only contribute to the de- 
polarized (VH) portion of the Raman spectrum while 
Eg modes only contribute to the polarized (VV) part. 
The Raman spectrum in our measurements is polarized. 
Thus, ujR is experimentally 3732 m/s and theoretically 
4837 m/s. 

A possible resolution of this poor agreement is that 
the gold nanoparticles are loosely attached to the silica 
matrix. In that case, FSM frequencies from Tab. |U would 
be appropriate. Since the E g mode frequency is 3779, 
the agreement with the experimental value of 3732 m/s 
is compelling. Furthermore, the next paragraphs show 
an independent reason to believe this. 

CFM calculations of the width of the peak show that 
it should be quite narrow. Based on the CFM value of 
ojR of the (SPH,2) mode given in Tab. [H] the expected 
FWHM width of a Raman peak at 22 cm -1 would be 
3.3 cm -1 , which is only slightly smaller than the ob- 
served FWHM width in Fig. ^ of 4 cm -1 . The increased 
peak width can be attributed to size variation of the gold 
nanoparticles within the sample as well as instrumental 
broadening. It is plausible to assume that the Raman 
width due to particle size variation, the Raman width due 
to mode damping, and the Raman instrumental width 
(1.4 cm -1 ) combine squarewise to give the width of the 
observed Raman spectrum. In that case, the Raman 
peak FWHM width due to particle size variation is just 
1.8 cm -1 , since 4 2 - 3.3 2 - 1.4 2 = 1.8 2 . We assume 
the distribution of particle radius is Gaussian with one 
standard deviation a and mean R aV g- The FWHM of a 
Gaussian distribution is 2V2 ln2 a. Thus, the rms varia- 
tion of particle radius is estimated to be a = 0.034R aV g. 

It is difficult to believe that the particle radius varia- 
tion in our sample is really so small, particularly since 
ion implantation experiments have never given such a 
tight size distribution. An alternative is that CFM over- 
estimates the true broadening due to damping. If the 
interface between the gold and the silica is not a tight 
bond then the broadening effect of the silica matrix is 
reduced based on a model of the interface as a thin soft 
layer^ Since gold does not readily form chemical bonds 
with silicon or oxygen the adhesion between the gold and 
silica is likely of Van der Waals type. 
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